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Abstract 

We discuss a two dimensional SU (N) Yang-Mills theory coupled 
to massive adjoint fermions for different worlds classified by the in- 
teger k = 0,1, ...N — 1. We study the fermion condensate for these 
unconnected worlds as a function of the parameter k. We show that 
the condensate as well as the spectrum of the theory do depend on 
this vacuum parameter k. 

Technically, the value of the fermion condensate is related to the 
value of the gluon condensate via the Operator Product Expansion. 
We use this to find the leading dependence (in the limit of a heavy 
quark) of the fermion condensate {k\ifjip\k) on the nontrivial vacuum 
angle k. We also determine the gluon condensate of the theory using 
low-energy theorems. 



1 Introduction 



Two dimensional SU(N) Yang Mills theories coupled to adjoint fermions 
continue to attract theorists 0-[ll| in view of their interesting connections 
with both string theory on one side and real four dimensional QCD on the 
other. In particular, as in the latter, QCD2 with adjoint matter is expected to 
possess a nonvanishing fermion condensate (ipip)- This corresponds to a non- 
trivial vacuum structure of the theory. It is also believed that such a structure 
is a consequence of nontrivial topological properties of the group SU(N)/Z N , 
which is the relevant symmetry group when adjoint matter (rather than fun- 
damental matter) is considered^. Therefore, it is expected that the theory 
possesses a nontrivial vacuum angle k |12| , similar to the 9 angle of QCD. 
It is known that many problems such as the U(l) problem, chiral symmetry 
breaking phenomenon, confinement and multiplicity of vacuum states are re- 
lated to each other and to 9 dependence. Therefore, in spite of the fact that 
we live in a certain vacuum state (this is expressed by the so-called supers- 
election rule), the physical parameters of the theory do depend on 9, which 
is an extremely important characteristic of the theory. 

The main goal of the present paper is an analysis of the condensate 
(k\%j)ip\k) as a function of the discrete topological vacuum angle label k = 
0, 1, ...N — 1. This dependence can be found explicitly as we shall see. There- 
fore, the vacuum condensate could be considered as an order parameter which 
distinguishes different vacuum states. 

First we would like to give a short historical introduction of the subject. 
For the case of fundamental fermions in the large N limit ('t Hooft model 
[ |l"3|j ) and vanishing fermion mass, the fermion condensate was first calculated 
in WM. Later on, the quark condensate was calculated in several different 



ways, both analytically and numerically QT5J [|T6| ||T7 



In the case of adjoint fermions, using bosonization it has been argued 
[0] that the condensate arises for arbitrary N. At the same time, instanton 
calculations seem to imply the vanishing of the condensate in the chiral 
limit for N > 3. An independent argument, based on quark-hadron duality, 
has been given to support a nonzero value for the condensate in the large 
N limit [El. On the formal side, in the small volume limit the condensate 



1 see however, Ref.jy] where it is argued that these topological properties are relevant 
for a theory with fundamental matter also. 
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was calculated for SU{2) and SU(3) gauge groups using ET quantization 
|J. Similar calculations have been performed for the SU(2) theory in finite 
volume L using light cone quantization (see 0). In all finite volume 
calculations the value of the condensate appears to scale as where L is the 
quantization volume. Therefore, no continuum limit can be taken in these 
formulae, and no finite result can be obtained. However, it is believed that 
an extra factor of L will appear in a complete theory leading to a finite value 
for the condensate as in the Schwinger model [Hj. Unfortunately, this 
would require a complete solution of the problem which is not possible at 
present and thus the question regarding the magnitude of the condensate 
in the continuum limit can not be answered within this approach at the 
moment. 

We conclude this short historical introduction by emphasizing that no 
meaningful calculation for the quark condensate is hitherto available in the 
continuum limit for two dimensional QCD with adjoint matter. Therefore, 
discussion of the number of different vacuum states which is based on the 
presumption of a nonzero value for the quark condensate has no solid basis 
before a finite value for the condensate (as an order parameter) is obtained. 

Before proceeding we should give a definition of the quark condensate in 
the theory. We define the fermion condensate in PCAC terms fll8l 0: 



This is the standard nonperturbative definition of the quark condensate in 
terms of the original tp fields. In particular, in the chiral limit m — > 0, only 

— £ 1 

a Goldstone boson with a finite (at m — > 0) residue (O|^>z7s'0(a;)|7r) = i ^ 
contributes to the correlation function (O|T'{'0z75^(j;)^z75^(O)|O). In this 
case eq.flID gives the famous relation m%f% = —4m(0\^'ip\0) between the 
mass of the Goldstone it meson and the chiral condensate. 

We want to emphasize that relation ([I]) is valid not only in the chiral 
limit m — > 0, but is also satisfied for arbitrary m. We expect, of course, that 
at large m — > oo the condensate goes to zero like (0|-?/^|0) ~ — . 

2 The singlet axial current in QCD2 is anomaly- free unlike the Schwinger model where 
there is an extra term due to the anomaly. 
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Nevertheless, a small, but not exactly vanishing condensate will play its 
role: it will be an order parameter which labels the different vac- 
uum states \k) we have been talking about. This is the cornerstone of our 
approach. We shall calculate the fermion condensate in the limit of large 
mass and we shall find different magnitudes for different \k). Presumably a 
calculation at large m corresponds to the weak coupling regime. Therefore, 
a perturbative calculation is justified. At the same time the small non- 
vanishing condensate is sensitive to the particular |fc)-vacuum state and has 
a different value for each vacuum state. 

As the first step in this program, we shall test our conjecture in the 
't Hooft model |13|]. This model is well understood in the limit of large N. 



We know the spectrum as well as all relevant matrix elements. Therefore, 
one can use definition ([!]) in order to calculate the condensate for arbitrary 
m exactly. The corresponding formula is also known and is given byJ17|]: 



mN 



I ren 

Z7T 



ln(7rn) - I - ~ L + ( 1 - - ) [( I -n)Z(n) - In i] 



(2) 



where the dimensionless parameter of the theory. We shall see 

in section 2 that a perturbative calculation is justified in the large m limitQ 
in this model. 

Encouraged by these results, we go one step further in section 3 and carry 
out the same perturbative calculations for the case of adjoint fermions with 
arbitrary iVQ. We next use the Operator Product Expansion in the limit 
m — > oo in order to relate the fermion and gluon condensates. Nontrivial, 
nonperturbative physics comes into the game through the gluon condensate. 
The dependence of the gluon condensate on the nontrivial vacuum structure 
has been determined earlier pj in pure Yang-Mills theory in 2 dimensions by 
means of the well developed machinery of Wilson loop calculus. 

Finally, in section 4, we use low energy theorems and the form of the 
fermion condensate we obtained to determine the -K corrections to the gluon 
condensate due to a large, but finite quark mass m. 



3 It corresponds to a small parameter a or small coupling constant g. 
4 Note that this model has not been solved yet, even in the large N limit. Therefore, no 
information similar to the 't Hooft case is available. 
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2 Lessons from the 't Hooft Model 



The 't Hooft model [13] consists of quarks interacting via gluons of the SU(N) 
gauge group with Lagrangian: 



£ = -\G; u G a ^ + $ (iYD, - m) V>, (3) 



where 



G% = d,A u a - d v A« - gr bc A, b A u c , (4) 
D^^d^ + igA^T^, (5) 

and T a are the Hermitian generators of the group representation. In the 
fundamental representation they are the matrices A a /2. This formulation of 
the theory is similar to fl9| and differs from the 't Hooft form in that here 
the gauge group is SU(N) instead of U(N). As mentioned in |L9] the U(N) 
singlet decouples and describes a free field and to leading order in 1/N the 
distinction does not matter. 
In light cone coordinates: 

x ± = ^=(x 1 ±x°), (6) 

the problem becomes very simple in the light cone gauge, A_ = 0. The 
Lagrangian becomes: 

£ = ^Tr(d^A + ) 2 + ij {yfDp - m) 4). (7) 

There is no ghost in this gauge. Take x + as the time coordinate and notice 
that A + is not a dynamical field but provides a (non-local) Coulomb force 
between the fermions. 

The algebra for the gamma matrices is: 

1- = 7+ = (8) 
7+7- + 7-7+ = 2. (9) 

The Feynman rules in this case are shown in Fig. ([J). 
Considering the limit N —>■ 00 with g 2 N fixed corresponds to only keeping 
the planar diagrams, 't Hooft solved for the dressed propagator in this case 
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' ( Y.P + + Y + P.+ m) 
U 2p + p - m 2 - ie 



a 





Figure 1: Feynman rules: where i,j are matrix indices of the group represen- 
tation and a,b are vector indices in the group space. 

using the Bethe-Salpeter equation. The form of the dressed propagator in 
our notation is: 



where A is an infrared cutoff parameter that removes the point fc_=0 from 
the momentum space, 't Hooft found that A drops out in all gauge invariant 
quantities. This regularization has been argued |19| |20| to be equivalent to 
a principal value regularization procedure where the photon propagator is 
replaced by it's Cauchy principal value (CPV): 



In this form A -1 is actually a gauge parameter which makes its cancellation 
from gauge invariant quantities obvious. With this regularization prescrip- 
tion the propagator is identical to the 't Hooft propagator except we neglect 
the terms involving A. 

With the motivation presented in the Introduction, we now want to cal- 
culate the value of the condensate in the well understood 't Hooft model in 
the large m limit. We achieve this by two independent approaches. First, 



z(7_p+ + + m) 



(10) 



2p + p_ — m 2 + g 2 N/n — g 2 iV|p_|/7rA — ie 



m 



Figure 2: Bare and dressed vacuum loop calculations shown to agree to order 
9 2 - 

we do the standard perturbative calculations which are perfectly justified in 
the large m limit. The second way of doing the same physics is to make use 
of the dressed 't Hooft propagator which is a solution of the Bethe-Salpeter 
equation. Agreement of the obtained results should be considered as a test 
of the method where a nonperturbative condensate (H) in the large m limit 
is simply determined by a pure perturbative calculation. 

The first approach is to calculate {i/ii/i) from the first diagram in Fig.(|2|). 
We calculate this bare diagram using both the 't Hooft regularization scheme 
and principal value regularization and found that they are equivalent at least 
to order g 2 \ 

1 g 2 _ (N 2 - 1) W 
it 2 m 2 IT 2 771 

where the second Casimir constant of the fundamental representation arises 
from TfjTji = Tr(C 2 (F)l F ) = C 2 (F)N after contraction over group matrix 
indices because of the closed fermion loop. We note that in the large iV limit 
where g 2 N ~ const., the condensate (ipip) ~ N as it should. We also note 
that in the large m limit it goes to zero as 1/m in agreement with the general 
discussions presented in the Introduction. 

The second approach to the same calculation is to take the dressed 
't Hooft propagator and expand it to the first order in g 2 (from the defi- 
nition of the condensate we have to take the difference between the dressed 



m = -N C 2 {F) - y - = -^—^L — — , (12) 
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and free propagator): 



dp-dp + i{l-P+ + 7+P- + m) 

(2tt) 2 [ U 2p + p_-m 2 + ^ 2 iV/7r-^ 2 iV|p_|/7rA-ie 1 j 
i(7-P+ + 7+P- + m) 



2p + p_ — m 2 — ze 



For small values of g we can expand the denominator of the first term so that 
we have: 



1 



2p + p_ — m 2 + g 2 N/n — g 2 N\p_ — ie 
1 



2p + p_ — m 2 — g 2 N\p_\/n\ — ie 



g 2 N/Tt 



2p + p_ — m 2 — g 2 N\p_\/it\ — ie 



The first term of the expansion cancels the bare propagator in CPV regular- 
ization where terms containing A are simply dropped. Therefore we have the 
integral: 

dP ~ dP+ Sij rn g2N/7 l (14) 



(27r) 2 [2p + p_ — m 2 — ie] 2 

The result is the same as above ( |T2"D in the large N limit and agrees 
with the large m expansion of formula (0), (see Ref . [jP7H ) . The moral of this 
section is that we are confident that the perturbative calculation gives the 
correct expression for the condensate ([!]) in this completely solvable example. 
Now we would like to apply this approach to a model which (unlike the 
't Hooft model) is not completely solvable, but is much more interesting 
because it includes nontrivial vacuum structure labeled by an integer number 
k = 0,1, ...N — 1. We consider this problem in the next section. 



3 Condensate with Adjoint Fermions 

Now we want to do the same calculation for two dimensional QCD where we 
take the fermions to be in the adjoint representation. The QCD lagrangian is 
the same as before with the only difference that the generators of the group 
are the structure constants of the group 

(T a ) bc = -if abc a = l,2...(iV 2 -l). (15) 
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We now repeat our previous calculation of the first diagram in Figure (0) 
with the following result: 

$y) = -(n 2 - 1) c 2 (A) ^4 = ~ N ( N2 ~ !) -r-. ( 16 ) 

where we take into account the value for the second Casimir constant for the 
adjoint representation: 

TgTfi = Tr(C 2 (A)l A ) = C 2 (A)(N 2 - 1). (17) 

This result corresponds to the trivial vacuum, i.e. k = 0. We would like to 
do the same perturbative calculations for a nontrivial vacuum state \k) as 
well. Recall that the case of fermions in the adjoint representation has long 



been known [T_2| to possess N different vacua corresponding to fundamental 
external charges at the boundaries of the universe. This is equivalent to 
a universe with a background color electric field analagous to the #-vacua 



in Coleman's analysis of the Schwinger model |^TJ. Different vacua in the 



nonabelian case have, in general, different vacuum energy and can be labeled 
by discrete values of 9: 

0k = ^k fc = 0,l,...(JV-l). (18) 

These vacua are unstable when the dynamical fermions are in the funda- 
mental representation because they can screen the charges at the edges of 



the universe to form a trivial vacuum \12\. This is the analogy of pair cre- 



ation in the massive Schwinger model. When the fermions are in the adjoint 
representation such a screening cannot occur and therefore there is a non- 
trivial vacuum structure and consequently a nontrivial dependence of the 
observables on this vacuum angle k. 

Now the question is: Is it possible to obtain a nontrivial depen- 
dence on the topological angle k from the perturbative calculation, 
which apparently does not contain any information regarding the topological 
properties of the theory at all? Are we able to see the effect of the vacuum 
angle 9 k in a perturbative expansion ~ J2n c n(^) n or is this impossible in 
principle? 

The answer to the question formulated above is: yes, we can extract in- 
formation about vacuum topological properties of the theory. In fact, this 
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information is contained in the coefficients c n of the perturbative expansion 
in powers of g 2 /m 2 . The point is that external charges at the boundaries of 
the universe cannot, of course, talk in the gauge theory to vacuum fermions 
directly, but they can communicate via gluons as soon as the external charges 
act as sources for the latter. It is therefore clear that in order to see a nontriv- 
ial k— dependence in the fermion condensate, {ipip), we should go beyond the 
leading term ([Tj]) in the expansion over g 2 /m 2 and take into account terms 
of higher orders in 1/m which are accompanied by powers of the gluon field. 
Technically, this problem amounts to a calculation of the heavy fermion loop 
in an external gluon field (see Fig.(|3|)) which does depend on the topological 
number k. In our case both the use of the perturbative fermion propaga- 
tor and retaining only a finite number of powers of the external field (and 
its derivatives) are justified by the smallness of the parameter g 2 /m 2 . This 
procedure results in a particular form of the Operator Product Expansion 
in powers of 1/m (known also as the heavy quark expansion). Note that by 
gauge invariance the expansion in external field starts from a term of second 
order in the field, and this is why we expect to get a nonzero effect at the 
level of order g A /m 3 . 

Therefore, at large m we are justified in the use of perturbation theory 
in order to reduce the problem of the calculation of (Jk\ipip\k) to the problem 
of the calculation of the gluon condensate (k\TrG IJLll G llv \k) . The correspond- 
ing method is well-known and we refer to a nice technical review |22|] of the 
subject. The result of the calculation is:0 

<£#V#> - (0|#|0) = - ((k\G%G a ^\k) - (0\G%G a no)) 

+ 0(l/m 5 ). (19) 

The problem of the dependence of the gluon condensate (k\G a ^ u G a,lv \k) on 
nontrivial vacuum structure in pure gluodynamics in two dimensions has 
been solved earlier in O. The dynamical heavy quarks with mass m will 
bring some order 1 / m 2 corrections to this result and we neglect these correc- 
tions at the moment (see the corresponding calculations in the next section). 
We quote the relevant formula for the gluon condensate in Minkowski space 

5 We note that the analogous formula in real four dimensional QCD4 takes the form (22): 
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Figure 3: Heavy fermion loop in an external gluon field 



rather than in Euclidean space, where it was originally obtained 0: 

(k\G%G a nk) -(k = 0\G%G a ^\k = 0) = g 2 *KN - k)(N + 1) ^ (2Q) 

Therefore, the final expression for the condensate as a function of the vacuum 
label k takes the following form: 

mm - (oi«-io) = N± M=m±£ . (21) 

07T mr N 

We close this section with a few remarks. First, as was expected, a non- 
trivial dependence on k appears in the final formula only at the level ~ g 4 /m 3 
and therefore is relatively small in large m limit. Nevertheless, formula 
demonstrates that all physical properties are different for different vacua la- 
bel k. Vacuum condensates are different. String tensions]?]] are different and 
therefore the mass spectrum is different. Vacuum energies are also different. 
We demonstrated this fact by keeping only the leading 1/m terms. However, 
we believe that the statement has a more general origin and it is hard to 
believe that different vacua could become identical for a finite value of m as 
long as the two parameters m and k are independent of each other. Our last 
remark regarding formula fliTS| ) is that the function which appears there has 
a property which was expected from the very beginning|12j - it is a periodic 
function such that states \k = 0) and \k = N) are equivalent. Therefore, 
we do not obtain a new state each time we increase the label k. Instead, we 
return to the starting state after N steps have been made. We notice also 
the symmetry property k (N — k) which was also expected p]]. 
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4 Low Energy Theorems 



The purpose of this section is twofold. First of all, the low energy theorems 
make it possible to investigate the behavior of vacuum condensates with 
changing quark mass m in the vicinity of parameters where those condensates 
are known, i.e. at m — > oo. This information might help in understanding 
the qualitative features of the model in general and its vacuum structure 
in particular. The second goal of this study is the calculation of the -\ 
corrections to the gluon condensate calculated in pure gluodynamics J7| due 
to the presence of matter fields. 

A similar study of the low energy theorems in the 't Hooft model was 



carried out for the first time in ref. [14|, for vanishing fermion mass, where 



the corresponding quark and gluon condensates have been calculated. This 



program was pushed forward in Ref. [17] where the previous result was gen 



eralised to arbitrary quark mass in the 't Hooft model. 

The main idea of the derivation of the low energy theorems is quite 
simple |p3fl . From the definition of the Path Integral, the variation of the 



gluon condensate with mass m is determined by the following correlation 
function 

A-(G%G a n = -ilimj d 2 xe^(0\T{^(x),G; u G a ^(0)}\0). (22) 

At the same time, the obtained correlation function can be explicitly calcu- 
lated due to the following identity^] 

ifirraj d 2 xe^(0\T{^(x),G%G^(0)}\0) = -l-A-(^). (23) 

Up to this point the derivation in QCD2 is identical to that in QCD4. The 
only relevant point is the form of the Lagrangian @ and not the specific 
properties of the theory which of course depend on the dimensionality of the 
space-time. The difference appears in the explicit calculation of the right 
hand side of Eq. fl2~3|) . 



In QCDi the corresponding calculation |23| of d r{j g i\ (V^) is based on the 



observation that due to the asymptotic freedom and to the absence of mass 



6 To derive this identity we simply rescale the gluon field — gA^ such that depen- 
dence on g appears in the Lagrangian only in the combination — ■^G'^ v G a >*" . 
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parameter other than M (the ultraviolet cut-off), the dependence on g comes 
exclusively through Aqcd ~ Mexp — b ^ M ^ where b = 11N/3 — 2JV//3 is the 
first term in the Gell-Mann-Low /^-function. 

In the present case of QCD 2 the dependence on g 2 is known exactly in 
large m limit (21). Therefore, we can explicitly calculate the behavior of the 



vacuum condensate {G a G a ^ lL ') with variation of the quark mass m: 



d (lra ra » v , , r a ™, \ N % 4 k(N - k) {N + 1) 

^{(G^G ) k -(G fll/ G )oj = 7 ^ ^ • (24) 

Note that in this formula we keep only the information relevant to the k- 
vacuum dependence. As mentioned above, there is a leading term (~ g 2 /m) 
in the fermion condensate (|1~6| ) which does not contain any k— dependence. 
A similar term in the gluon condensate depends on subtraction and nor- 
malization procedures but we will not discuss this issue here. However, the 
difference in condensates between two vacuum states has an absolute mean- 
ing and the corresponding expression is given by eq. (|24|). 

Formula (|24]) makes it possible to calculate the finite correction due to 
the dynamical heavy quark to the gluon condensate calculated previously [0 
in pure gluodynamics. Integrating Eq.(|24D at large m, we obtain 



(k\G%G^\k) - (0\G^\0) = 9 2 2HN l )(iV + 1) (l - \ 9 — 2 )- (25) 



The correction ~^in this formula is not suppressed in the large N limit. 
This is in accordance with the general expectation that fluctuations of adjoint 
matter fields (unlike fundamental matter) are not suppressed by a factor 1/N. 

One can go a little bit further in the analysis of the low energy theorems by 
taking advantage of the specific property of two dimensional QCD that the 
only dimensionless combination in the theory is g 2 /m 2 . Therefore, without 
loss of generality one can present the fermion condensate in the following 
form: 2 

(^)=m/(-^), (26) 



.,2 



with some function f(x = -^). In this case the derivative with respect to m 
and derivative with respect to g are related to each other. At the same time, 
as we observed earlier, the derivative d /y g i\ (^V^) is reduced to some correla- 
tion function (|23"D which describes the dependence of the gluon condensate 
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on m (p2). Therefore one can relate the gluon and quark condensate exactly 
without any approximations using only very general properties of the theory. 
Indeed, from identities 




m° 



■/'(*), 



one obtains: 



d(l/gi 



g m d ( (ipip) 



2 dm 



m 



(27) 



(28) 



Now combining the low energy theorems 



exact relation between quark and gluon condensates: 



2 dm p dm \ m 



2|) with (|28|) we obtain the 

(29) 



This relation was derived earlier in |L7j using a different approach. One 
can rederive our previous formula for the gluon condensate ([25]) using the 
exact relation ([29]) and the expression (pl|) for the quark condensate in the 
asymptotic region. 

An exact relation between quark and gluon condensates should not be 
considered as a big surprise. Such a relation is a direct consequence of the 
fact that in two dimensional QCD the gluons are not dynamical degrees of 
freedom. 



5 Conclusion 

The main results of this paper are given by Eqs.(PT|),(P5]). We explicitly 
calculated the quark and gluon condensates in QCD 2 (N) coupled to adjoint 
matter as a function of the nontrivial vacuum label k. Our formulae are 
valid only in the vicinity of large m, but the obtained results suggest that all 
observables do depend on the specific vacuum state in which we live. 
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In principle, one can further generalize formula (|19|) in order to include 
higher dimensional gluon condensates. One can derive the next term in the 
expansion: 

m = -^r-Tr [T(R)T(R)] - ±-^(Tr [CTG^]) 

<r f (Tr [(D a Gn(D a G^)]) + \{Tr \G. V D 2 G^\) (30) 



127rm 5 V LV ° yv 2 
- \{Tr\{D a G^){D^)\ )) + ... , 



where T(R) are the generators in the representation R of SU(N). Note that 
Eq.(]30f) is valid for arbitrary representation of SU(N). Unfortunately, we do 
not know the higher dimensional condensates which appear in (|30| ) . The only 
thing we know for sure is the fact that the fc-dependence of each condensate 
should be a periodic function similar to, but not necessarily the same as, 
fl2ip. One can reverse this argument. If we knew the fermion condensate in 
QC D^ mnt for arbitrary m exactly, (similar to Eq.(||]) in the 't Hooft model), 
we would be able to calculate all high dimensional gluon condensates as a 
function of k. This problem is probably too complicated and equivalent to 
the complete solution of QCD^ mnt which presently is not available. 

We conclude this paper with the following project for future investiga- 
tion. One may try to use our results with heavy quark mass to generalize the 
string picture originally derived in |24] for the pure gauge 2D theory without 



dynamical degrees of freedom. By introducing a heavy quark we have essen- 
tially introduced the physical degrees of freedom without noticeably changing 
the internal gluodynamics. Therefore, we expect that the introduction of the 
heavy quark might be the first step in the direction of a string description 
of dynamical degrees of freedom. We probably need to implement some dy- 
namics on the string-sheet boundary which would correspond to a heavy, but 
dynamical, quark. 
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